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f—>. ' We study the level of primordial non-Gaussianity in slow-roll two-field inflation. Using an analytic 

f— ») , formula for the nonlinear parameter /nl in the case of a sum or product separable potential, we 

pvj ■ find that it is possible to generate significant non-Gaussianity even during slow-roll inflation with 
Gaussian perturbations at Hubble exit. In this paper we give the general conditions to obtain large 

£_) , non-Gaussianity and calculate the level of fine-tuning required to obtain this. We present explicit 

f~\ . models in which the non-Gaussianity at the end of inflation can exceed the current observational 

^ ■ bound of |/ NL | < 100. 

^O . 

PACS numbers: 98.80.Cq 

E" 
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£ ■ I. INTRODUCTION 

GC 

Standard slow-roll single-field inflation generates a quasi-scale invariant spectrum of Gaussian, adiabatic perturba- 
tions. While this is in agreement with current data, string theory and SUSY gencrically contain many scalar fields 
and it is important to test how these extra fields may change the predictions of the simplest models. In particular, 
the simplest models predict a level of non-Gaussianity which is too small to observe in the foreseeable future [l|, 0, S] , 
so any detection would be of great significance. 

In this paper we focus on the possibility to obtain a large level of non-Gaussianity during slow-roll inflation and 
derive general conditions for when it may be large. We present simple, explicit two-field models which can saturate 
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the observational bound. Although there are many models of inflation which generate a large non-Gaussianity [4 
a, [D, |8|, |9|, [lCJ, few of them do so during inflation except by breaking slow-roll, e.g. with a kink in the potential ^ 
or having a non-standard kinetic term, e.g. DBI inflation [13] • Other widely considered methods to generate a large 
non-Gaussianity include the curvaton scenario [l3|, modulated reheating [141 ] or an inhomogeneous end to inflation 
fla Ha . fl7l [la . [l9| . All these scenarios require an additional light scalar field that doesn't affect the dynamics during 
inflation but that becomes important either at the end of inflation or later. 

Several authors have considered the possibility of generating non-Gaussianity during inflation for a separable po- 
tential. In particular Vernizzi and Wands calculated a general formula for the non linearity parameter f^L which 
parameterises the bispectrum in the case of a sum separable potential [2fJ, this was later extended to the case of a 
product separable potential including non-canonical kinetic terms [2l| and a sum separable potential for an arbitrary 
number of fields J22JI . This has been further extended to the trispectrum (4-point function) [23j]. For other approaches, 
see for example [2J, [25( ■ 

Rather less work, however, has been done on analysing the formulas which have been calculated and are present 
in the literature. The results are rather long and depend on many parameters. In general they appear to give a 
slow-roll suppressed non-Gaussianity, subject to several simplifying assumptions that are additional to the slow-roll 
assumption and that are not always valid. Here we carefully consider several explicit models of two-field inflation 
and scan the parameter space more generally than has been done before. We show that it is possible to generate 
an extremely large non-Gaussianity during slow-roll inflation even when the field perturbations at Hubble exit are 
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Gaussian. However this is only possible for specific values of the model parameters and initial conditions. We consider 
under which conditions a general model can generate a large non-Gaussianity and show that this always requires some 
fine tuning of initial conditions. 

Although previous papers have shown that it is possible to generate a narrow "spike" of large non-Gaussia nity 
while the inflaton trajectory turns a corner, the non-Gaussianity decays again very quickly after the corner [2fJ, |2J]. 
This spike of non-Gaussianity is associated with a temporary "jump" in the slow roll parameters [251 ]. In the models 
we consider the non-Gaussianity is large over many e-foldings of inflation and all of the slow roll parameters remain 
smaller than unity. 

The paper is organised as follows: First we define relevant quantities and introduce some notation (Section|TT|. Then, 
in the following two sections we derive the general conditions to generate a large non-Gaussianity during inflation, in 
Section IIIII for a product separable potential and in Section IIVI for a sum separable potential. In Section [V] we give 
specific examples for product and sum potentials which can generate large non-Gaussianity. This includes a two-field 
model of hybrid inflation in Section IV CI In Section IVII we extend the previous results to a generalised action with 
non-canonical kinetic terms. We conclude in Section lYlII 

II. BACKGROUND THEORY 

In this paper, we consider an inflationary epoch driven by two scalar fields, whose dynamics are governed by the 
action 
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Here Mp = 1/\/8ttG is the reduced Planck mass. We consider slow-roll inflation, during which all of the slow-roll 
parameters defined below arc less than unity. 
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While it is physically interesting to consider the slow-roll parameters, it will be useful to use the definition of angles 
along or perpendicular to the background trajectory of the two inflationary fields [2q : 
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such that, in the slow roll approximation, 
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Throughout this paper we use formulae given in [2l| and we refer the reader to this paper and references therein for 
more details. In summary we define several observable quantities in terms of the primordial curvature perturbation 
£, which may be calculated using the <5N formalism [271, S, S, 120, [3JJ . The number of e-foldings, N, given by 



N 



H{t)dt, 
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is evaluated from an initial flat hypcrsurfacc to a final uniform density hypersurfacc. The perturbation in the number 
of e-foldings, SN, is the difference between the curvature perturbations on the initial and final hypcrsurfaccs. In this 
paper we take the initial time to be Hubble exit during inflation, denoted by t*. In this paper, we consider the final 
time to be a uniform density hypersurface just before the end of inflation, as in [2fJ,l2]J. We note that this differs from 
the hybrid scenario of [18|. in which the final hypersurface is determined by a separate waterfall field. The curvature 
perturbation is given by [31( 
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where N,i = dN/(d(fl) and the index / runs over all of the fields. We will consider the power spectrum and bispectrum 
defined (in Fourier space) by 

2tt 2 
(C kl Ck 2 ) = (2^) 3 <5 3 (k 1 +k 2 )-^P c (fc 1 ), (7) 

(C kl Ck 2 Ck3> = (27rfS 3 (k 1 +k 2 +k 3 )B c (k 1 ,k 2 ,k 3 ). (8) 

From this we can define three quantities of key observational interest, respectively the spectral index, the tensor- 
to-scalar ratio and the non-linearity parameter 
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Here T 5 * is the power spectrum of the scalar field fluctuations and Vt = SV* = 8H% / (4ir 2 Mp) is the power spectrum 
of the tensor fluctuations. Both the spectra are calculated at the end of inflation and we ignore any evolution after 
this time. For full evolution after inflation, see [2l|, [321]. As defined above, /nl is shape dependent, but it has been 
shown that the shape dependent part is much less than unity [20, [33| • The ideal CMB experiment is only expected to 
reach a precision of /nl around unity |34j |. so we will calculate the shape independent part of /nl, denoted by / NL ^ in 

[20, [21(. Whenever the non-Gaussianity is large, |/nl| > 1, we can associate / N l — /nl- This (k independent) part 
of /nl and the spectral index can be calculated by the SN formalism, 
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There is no universal agreement over the sign of /nl in the literature. We use the opposite sign convention to [20, l2ll ] , 
in order to be in agreement with the WMAP sign convention [35[ . The latest observations from 5 years of WMAP 
data are 

n c = 0.96to.oi3, (assuming r = 0), (15) 

r < 0.2 (95% CL), (16) 

-9 < /j^ 1 < HI (95% CL). (17) 

The "local" bound on /nl is based on the definition C = Cg + 3/nlCg/5, where £g is the Gaussian part of £ |34.l3a|. 
When more than one field direction during inflation contributes to £ in Eqn. |6]) then this definition of /nl is not 
equivalent to the very commonly used definition Eqn. pip , even though they have the same shape dependence (see 
e.g. [36| )■ The bound on f^ al may therefore not be applicable to the /nl that we calculate. In this paper we will 
focus on the non-linearity parameter /nl- Constraints on this parameter are expected to improve by nearly an order 
of magnitude with the PLANCK satellite [37|. 

III. LARGE NON-GAUSSIANITY IN PRODUCT POTENTIALS, W(tp,x) = U(ip)V( X ) 

We follow the terminology of [20] an d use the definitions, 

u = ^ = cos 2 9 e , v=^= sin 2 9 e . (18) 



We will use a sub- or superscript "*" to denote values evaluated at the time of Hubble exit and a subscripted "e" will 
denote the time at t e . From Eqns. ()12p and (|13[) . the power spectrum and spectral index for a product potential with 
canonical kinetic terms are 
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The non-linear parameter f^J becomes [21 
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The slow-roll parameters, defined by Eq. ([2]), are 
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Using a few auxiliary functions (involving only 0* and e ), we may write Eqn. (|21[) as 

/nl - jj [2j P (r , e e )£* - uw^Ww - 9 P (e*,e e ) v * xx + 2h p (e*,e e , x) (f - 

where the auxiliary functions are defined as 
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By analysing the functions j p , f p , g p and h p over the range of allowed values for 0* and e , it is possible to locate 

regions of parameter space which gives large / NI \ We find that the function j p (0*,9 e ) satisfies < j p < 1 and 

therefore the 'j '-term can never lead to significant values of / NI ^ ■ We therefore ignore this term in the analysis that 
follows. 

The functions f p , g p and h p (as functions of 0* and e ) are plotted in Fig. [TJ The main point to note is that these 
prcfactors are large (and can be very large) in two regions: 
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FIG. 1: The contour plot of the functions, f p , g p and h p , in the plane of 9* , and 9 e . The bottom and left-hand axes show the 
angles, 9* and 9 e respectively. The top and right-hand axes show sin 2 9* and sin 2 9 e . 



A cos 9* <C 1 and cos 9 e <C 1 (e x ^> e^,). In this region, f p ^$> g p and h p ~ / p sin 6* e 
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In both these regions, these prefactors can account for a significant level of f^, despite the relevant terms in 
Eqn. (|26p having slow-roll factors. Each of these regions describes one of the fields dominating over the other in 
kinetic energy. Note, however, that single-field inflation (i.e. with one field static) will not lead to a large value of 
non-Gaussianity, since these fractions vanish exactly for x e — or ip e — (6 e = or w /2 respectively). 

Due to the symmetry, we shall focus on Region B and explicitly write down full conditions for both Regions A and 
B in Section IIII Al To find the condition to give large non-Gaussianity, we will concentrate on the terms including 
the functions of g p (0*, 9 e ) and h p (9* , 9 e ), since the other terms cannot give large non-Gaussianity as discussed above. 

The large f^l is given by 
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Here we assumed that \r] xx \ ^> sin 2 9 e \rf vv — 4e e |. For |/nl| to be larger than unity, the function g p (9*,9 e ) must be 
bigger than the inverse of the slow- roll parameters in the square bracket. In this way from Eqn. (|28[) we find the 
general way to obtain large non-Gaussianity, \f^\ -^ 1, 
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This bound is shown graphically in Figure [2j 

From this condition we can derive some corollaries. Firstly, from dg p /d(sin 2 9 e ) — and the definition of g p , the 
initial and final angle of the trajectory are constrained: 
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and secondly we see that sin 9 must typically grow by at least two orders of magnitude during inflation, since we 




FIG. 2: A blown-up graph of Region B. The conditions for g p in Eqn. (|29[) is plotted for Q p = 1000 with a white line. It can 
be seen that this condition encloses the contour for g p larger than some constant, Q p . 
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In this region of large non-Gaussianity, the power spectrum from Eqn. (|19[) . the spectral index from Eqn. (|20p and 
the tensor-to-scalar ratio in Region B are: 
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The spectrum of curvature perturbations can be dominated by the fluctuations in the slow-rolling field ip or \ 
depending on f, which can be smaller or bigger than unity. We note, however, that the tensor-to-scalar ratio can only 
be suppressed in multiple field inflation, since the tensor power spectrum is unchanged from the single field value and 
Vq is enhanced. The spectral index is changed by the presence of f , but remains slow roll suppressed. 



A. Condition for Large Non-Gaussianity in Product Potentials 



If we summarise, the general conditions to obtain large f-^J can be written simply. Ultimately, the two-field system 
must have one field which dominates the evolution throughout. This corresponds to being in Region A or B (depending 
on which field dominates). Hence, one field must dominate the evolution almost fully at the beginning, while the other 
field gains a relatively large percentage of kinetic energy by the end of inflation. Specifically, the slow-roll parameters 
and the value of f p or g p must be such to give |/nl I ~ 1 m Eqn. (f2l)|) . 



For Region A (cos 2 9 <C 1), the condition is 
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For Region B (sin 2 9<1), the condition is 
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This condition also leads to a large growth of the subdominant slow-roll parameter. For Region A and B respectively: 
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B. Direct Observation of /^ and Fine-Tuning 

Eqns. (|36|) and (|37|) encode the level of fine-tuning required on the initial and final slow-roll parameters to obtain 
an observable level of non-Gaussianity. As an example, if observations find fLJ ~ 10, for standard order of slow-roll 
parameters (77 ~ 0.01), one requires f p or g p ~ 1000. In the case of Region B, this corresponds to 

A = sin 2 8* < 10~ 7 , and A = sin 2 6> e < 10~ 3 . (39) 

This requires very special initial values for the fields. 

IV. LARGE NON-GAUSSIANITY IN SUM POTENTIALS, W(ip,x) = U{f) + V(x) 

In this section we find the general condition for large non-Gaussianity with a sum separable potential, W(ip, x) = 
U(<p) + V(x). Defining 
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the power spectrum and spectral index are given by [20j: 
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The nonlinear parameter fLJ is [201 ] : 
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The slow- roll parameters, defined by Eqn. ([2]), are 
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Similar to the analysis of a product potential, it is possible to re- write Eqn. (|44|l : 
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where the following auxiliary functions have been used: 
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Note that u + v = 1 and hence these functions depend on just two variables, u and 9* (or v and 0*), and that 
< u < 1 and < u < 1. The first term of / NI ^ m Eqn. (|4"9")l is always smaller than unity, since < j s < 1. 
Once again there are two regions with potentially large non-Gaussianity: 

Cu«l, and cos 2 6* < 1 

D v < 1, and sin 2 0* < 1. 

We temporarily analyse this second region, Region D, in which v <C ljSin 2 ^* <C 1 and g s can be large. In this 
region, we can approximate 
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In the special type of potential when V e ~ V*, such as hybrid inflation, v ~ sin 9 e W e /W* and hence 
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The function g s has exactly the same form as the function g p for the product potential in the region where it is 



large, see Eqn. (|28|) replacing sin 2 8 e with v. Hence it is large in exactly the same areas, and for g s > Q s we require 
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Finally, we can approximate /^ in Region D as 
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where we have used u ~ 1 and cos 2 8* ~ 1. We may analyse this equation for the specific case when V e ~ V*, which 
can be true in many models because sin 2 8* ~ (V x /U v ) 2 <C 1, so at least initially the potential V is extremely flat. In 
this scenario, v ~ sin 2 8 e W e /W* and hence 
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Indeed, the form of power spectrum and spectral index in this large non-Gaussianity region, are same as that of 
product potential, Eqn. (f3"2"|) and (|3U|) . 

A. Condition for Large Non-Gaussianity in Sum Potentials 

From the previous section, it is clear that, to get large non-Gaussianity, we need a constraint on the function g Sl 
for Region C, 
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In the same way, for Region D, we obtain 
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As in Eqn. (|30[) , wc find exact upper limits on the parameters: 
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V. EXAMPLES 

In this section we give specific examples, two of which can give large non-Gaussianity and one which cannot. We 
note that, in all the following examples, we have one parameter that acts as a normalisation factor in W 7 *. This 
parameter can be fixed to normalise the amplitude of the perturbations. 



A. Quadratic times exponential potential, W(ip, x) = 2 e 
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We consider a product potential, with U(<p) — e Xip l Mp and V(x) — m 2 x 2 /2, for which the slow-roll parameters 
are: 
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FIG. 3: The contour plot of log 10 |/^| for Example A, W(<p,x) = \e~ x ^ IM rrnx ■ Here we used A = 0.04 (left) and 0.05 
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(right). White regions indicate when \5^{\ < 1 and the dark centre-most region indicates I/^lI > 100 
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We consider A > and slow-roll throughout inflation, so that the yj-field increases during inflation while the x-field 
decreases. In this slow-roll limit, the exponential function can be expanded so that the potential is dominated by the 
quadratic potential of \ field, ■km 2 X^_ l and corrected by the integration between \ an d ¥> field, — ^Xm 2 tp 2 x 2 , which 
might be important for preheating |38l.l39j. Since inflation ends when the \ field rolls close to its minimum, i.e. when 
e = 1 ~ e x , the only way to generate large non-Gaussianity is to start with small field of ip> (which corresponds to 
Region A). 

For this potential, the slow-roll solutions for tp and x lead to 



Since X 2 (p 2 /Mp < 1, we find rf v v 
can be obtained: 



„2AJV 



ANMp. 



= <p*e— x~ =x;-4i\M P . (61) 

-2 A. From Eqn. (]36[) . the constraints on the initial values of the field tp 



where N e ~ (xl - Xe)/ 4M J> an d X 
constraints lead to 



2,.4_-12AJV e 



M P xte 



< ^< 
~ Ml 



IMle 



4AAT 



(62) 



2A3 X f ~ M 2 p ~ Xxl 

2 ~ 2. For the specific values of x* = 16A/p and A = {0.03,0.04,0.05}, the 



{0.134, 2 x 10~ 3 ,3x 10" 



'} ^ inr ~ {0.128,0.031,0.0078}. 

Alp 



In this range, the non linearity parameter is well approximated as 



f (4) 
'NL 



COS 6 e e 



6 (cos 2 9* + cos 4 



'Vw 



(63) 



/ A\ _____ 

The contour plots of J^l f° r ^ w0 values of A (A = 0.04,0.05) arc given in Fig. fusing the full formula Eqn. (|44[) . 
which is almost same as Eqn. (|63p . Here it can be seen that /jU can be very large when ip* ~ 10 _3 A/p (A = 0.04) 
or tp* rv 10~ 4 A/p (A = 0.05). We note that the analytic constraints differ slightly from the plots, due to the value of 
rj vlp used. In the analytics, we assume r\ lpv ~ — 2A, but use the exact form (Eqn. (fST))) ') in the plots. 

The time evolution of /^ from the 8N formula is given in Fig. [51 where the final point is identified as some time 

during inflation. This method is also used in Vcrnizzi and Wands [2(| for analytic evolution. We can see that I/^lI 
increases sharply around N = 30. This can be understood from the evolution of two fields: from this time, the velocity 
of tp (or cos 2 9) increases (as seen in the third and fourth plot of Fig. [4j to satisfy the condition Eqn. ((36)) to obtain 

large non-Gaussianity. Then /^ begins to decrease around N — 45 as cos 4 9 e becomes bigger than cos 2 9*. Note that, 
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FIG. 4: The analytic evolution of /^ L and e x and c v for Example A. The evolution of both fields, ip and x are a l so shown. We 
used A = 0.05 and ip* = 10~ Mp, x* = 16Mp. We numerically solved the full equations of motion until e = 1. 



around the region for large |/nl|, the curvature perturbation of tp-ficld becomes comparable to, and then dominates, 

that of the x-hcld. Around A = 58, the x-field begins to roll so fast that cos 2 9 decreases and |/ N l| increases once 
more. However during this fast-roll stage, slow-roll breaks down and our analytic formula is not valid any more. Thus 
we calculate final values when e = 1. 

We note, for reference, that in this regime of large non-Gaussianity, the background trajectory must be at almost 
entirely in the one-field direction. We therefore expect reheating to take place in almost the same way as it would in 
normal chaotic inflation. 



B. Quadratic product potential 

In this section we show that a direct coupling between tp and Xi leaving to the product potential 



(64) 



cannot generate a large non-Gaussianity during inflation. From the equation A = J Hdt it follows that the fields 
evolve according to 



£ 



yl = yl- 4NMJ 
and similarly for \. Hence we have in Region B, where sin 2 9* <C 1 and 4NMp < <p\ <C x 2 , 



,2 ae 



T% 



2 - 4A^M 2 



6 e xl - 4AM 2 " xl 



P <^ =sin 2a* 



(65) 



(66) 



Hence sin 2 9 decreases during inflation and this model cannot generate |/nl| > 1- 



C. Hybrid inflation 



We consider a model of 2 field hybrid inflation. 

W{fp, X ) = 



Wn 



1 + a TTT + fllTPT 

Ml M 2 



(67) 



which is vacuum dominated, i.e. which satisfies jav? 2 ] ^C M^, and \Px 2 \ ^ Mp. Here we assume that inflation ends 
abruptly by another waterfall field which we don't write down in the potential above. We note that /nl may change 
depending on the details of how the waterfall field is coupled to the two inflaton fields [l5|, [l6|, [l~8[ , but in general this 
is unlikely to generate a large contribution to /nl without fine-tuning [17] . 
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0.018 -0.018 
0.04 0.005 
0.01 -0.02 



p*/Mp x*/M F 



0.00018 
0.0018 
0.00037 



./] 



NL 



-42 
9.27 
-11.1 



n^ — 1 




0.09 
-0.02 



0.006 
0.10 
0.026 



TABLE I: Table showing some initial conditions for the hybrid inflation model that lead to large levels of non-Gaussianity. The 
spectral index, calculated from Eqn. (|43p . and the tensor-to-scalar ratio are also shown. 

In this regime the slow-roll solutions are [5J] , (we can identify rj v = 2a and r\ a = 2/3 in the notation of that paper) 



if = tp*e 



-2aN 



X = X*e 



-2/3 JV 



and the slow-roll parameters arc 



?] vv = 2a, r] xx = 2/3, r) vx = 0, 

2 2 



(68) 



(69) 



We note that the dominant slow-roll parameters 77^ and r\ xx are constant in this model. 

For this case, since W e ~ W*, from Eqn. ([54")) (or equivalently Eqn. (jBT))) ) in Region D, where sin 2 9* <C 1 and 
v ~ sin 2 9 e < 1, it follows that 



/^^75,[-%x-2( e e -r)] = 



Sin 



6^ L ,AA ' /J 6(sin 2 6»*+sin 4 6l e ) 27?xx ' 

In the above we have used Eqn. (|51|) for g s in Region D. We hence require the condition, Eqn. (|57p . 



sin 2 (9* < u 2 



1 



- 1 



Gs = -^\Vxx\ \ 



and find the exact upper limits on the parameters: 



Sin20 * < l\l) (4) l%x|Z and v = sin ' 0C ■ ^Vxxl 



-■-^ r- ^ 



(70) 



(71) 



(72) 



Noting that in Region D sin# = ftx/( ai f)j from Eqn. (|6"8"]) we require N(a — /3) > 1 so that sin grows significantly 
during inflation. 

In Tabic HI we give some explicit examples of values of a, j3, (p* and x* which lead to a large non-Gaussianity. Using 
Eqn. (|4"3")) we also calculate the spectral index. The contours of /^ of this sum potential for a specific choice of a and 
P is given in Fig. O The first example in the Table [T] shows that it is possible to have /nl — 50 and a scale invariant 
spectrum. The tensor-to-scalar ratio is also small in this example. 

Our formula for f^£ is also consistent with the calculation of 0] but our analysis of the parameter space is more 
extensive and in particular it allows for very small initial field values. We therefore find a region with much larger 
non-Gaussianity than was found in Q ■ 

As a further check on the algebra, we note that it is also possible to analyse this model using the formalism 
of a product potential with W(ip, x) = Wo exp(avj 2 /M|>) exp(/3x 2 /M|,). This is equivalent to the sum potential 
Eqn. (|6T|) in the limit of vacuum domination (\a\ip 2 /Mp -C 1 and \(3\x 2 /Mp •C 1). We have checked that this 
gives the same results. We note from Table Q] that we require a very small value of x* m Region D. If the value 
becomes too small then the motion of the x field w iU become dominated by quantum fluctuations rather than the 
classical drift down the potential, 3Hx — —W x = —WqPxi which we have assumed. In order that we can neglect the 



effect of the quantum fluctuations, the condition we require on the background trajectory is \x\tt/H > \/3/2 [4 



Using Vc, — 10 10 , H 2 ~ Wo and Eqn. (|42[) the condition that the classical trajectory is valid can be rewritten as 
(1 + r) P 2 xl / (ct 2 ipl) > 6V(, where r = sin e /sin 2 9*. We have checked that this is satisfied for all of the examples 
given in the Table [H but it does provide a significant constraint on the total parameter space. 

Very recently Cogollo et al. have calculated the effect of the loop correction to the primordial power spectrum and 
bispectrum [411 ]. This loop correction arises from taking into account the contribution to the power spectrum and 
bispectrum arising from terms in the 6N expansion which are non-leading in the expansion of the field perturbation 
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FIG. 5: The contour plot of log 10 |/^lI with the sum potential, W((p,x) 

match the examples in Tableland the contours for log 10 I/^lI are shown. White regions indicate when I/^lI 



Vb(l + aip 2 + flx 2 ) (Example C). The parameters 

< 1. 



Sip, see e.g. [42|]. However they can still be significant if the coefficient to the term given by the derivative of N is 
extremely large, e.g. J43|. These "higher order" terms are usually neglected, but [41| has shown the first explicit 
example of an inflation model where they cannot be neglected. They consider a 2-field hybrid model with the same 
potential as Eqn. (|67p in the special case of an unstable straight trajectory along one of the axes, with a and (3 
both negative. In this case, they find (for certain initial values) the loop correction is dominant and can generate an 
observable /nl- We have checked that for the explicit values given in the table the loop correction does not dominate 
(under the assumption that the same loop correction is dominant here as the one in [4l|). However this does provide 
a further restriction on the allowed parameter space of the "tree level" calculation of the model, which is assumed in 
Eqn. (|70p . We plan to return to this issue, and make a more thorough investigation of the hybrid model in a future 
publication. 



VI. NON-GAUSSIANITY WITH NON-CANONICAL KINETIC TERMS (PRODUCT POTENTIALS) 



Generally, the scalar fields need not only couple through the potential, but may also couple kinetically [44l . l45l .l46l.l4'i 
In this section, therefore, we consider inflation governed by the following generalised action: 



5* = 



Mi 



R 



- Ig^d^ip - l -e 2b ^g^d^xd„X ~ W(<p, X ) 



(73) 



When the kinetic energies are non-canonical, the local non-Gaussianity is altered from the previous sections. In 
the case of the sum potential, the modifications are non-trivial and are difficult to analyse. When the potential is 
of product form, however, we may at least partially use the previous analysis. We therefore concentrate only on a 
product potential, but our findings relate also to sum potentials. 

With a product potential, the local non-Gaussianity is given by |2l| : 



6 , (4) . 
c-'NL — 



2e -2b c +2b, 



C* 



1 _ ? w 
2e* 



+ 2 sign (6^) sign I -^ 



(C) 2 



et \ 2e* 



tin 





26.-26, 



At 



(74) 



where 



n 



-2b e +2b, 



1+^(1- 



2&„-2b» 



^ps- 



*[l + tan 2 e (l-e X )], 



e e 

V> X 



I], 



^sign(6 v )sign ( -£- » 



[/ 
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(75) 
(76) 
(77) 



Using new auxiliary functions, Eqn. (|74|) can be re-written as 



/$ = [2J(r,0 e ,X)e*-F(r,0 e ,X)r^-G(r,r,X)ry* x + sign (&„) sign (^JK(0*, e , X)^*^ 



-2H(6*,e*,X) rf- isign(6 v )sign^J ^- 



4 /,e 



„ J e e -i/^-4sin 2 e COS 2 0V 



(78) 



where 



F„ 



(~ T i. 



H„ 



.X) 

~x) 



~ x a 3 tan 4 0* 



(a 2 tan 2 0* + tan 4 6» e ) cos 2 e 

e~ A " tan 8 e 
(a 2 tan 2 0* + tan 4 e f sin 2 6> e ' 

,2 /)* + „„2 /)e\ 2 



\0 e ,X) = tan 



2/ae 



(a tan^ 0* - tan 2 
(a 2 tan 2 0* + tan 4 6» £ 



>2 • 



j p (r,e e ,x) 



K p (0*,0 e ,X) = e 



-X 1 



0* (a 3 tan 2 0* + tan 6 e ) 
cos 2 e ( Q 2 tan 2 0* + tan 4 e V 



-X„2 



cos 4 e sin 2 i 



sm 



(79) 



If X = 0, then a — 1 and i^ = / p , G p = g p , if p = /i p and J p = j p . The function F p is plotted for various values of 
X in Figs. [5] From this definition, it is clear that allowing a =/= 0, the range of e for which we can obtain large f NL 
opens up. We also note that the symmetry between two fields is broken, due to b(ip), which is apparent in Figs. [5] 
Note, though, that we require 0* ~ 0, ^ as before. Furthermore, when X > (or b e > b*), we no longer require 
cos 2 e < 1. 

Expressions for the power spectra and spectral index with a non canonical kinetic term are given in [2l|, [47] . 



VII. CONCLUSION 



We have made an in depth investigation into the level of non-Gaussianity during two-field slow-roll inflation. We 
have shown that it is possible to generate a large level of non-Gaussianity during inflation without violating slow roll 
and when the inflaton field perturbations arc Gaussian at Hubble exit. 

The general conditions for generating a large non-Gaussianity show that the inflaton potential must have a specific 
shape so that the angle of the background trajectory can grow by about two orders of magnitude or more during 
inflation. In the case of an inflation potential made of a product of two quadratic potentials, this condition is not 
possible so that we conclude this model cannot generate a large non-Gaussianity during slow roll inflation. When the 
angle of the background trajectory grows sufficiently (in relative terms), we still need one of the fields to dominate 
throughout inflation, yet the remaining field can not remain full insignificant. For typical values of the slow roll 
parameters, we require that initially the sub dominant field, say ip, must satisfy cos 2 0* — e*/e* < 10 . This means 
that the field trajectory is almost exactly parallel to the x axes initially and typically this requires a finely tuned 
initial condition. In the case of a product separable potential we then also require that the final value of the angle of 
the background trajectory lies in a narrow range much greater than the initial value (in relative terms) but still nearly 
parallel to the \ axis. The analysis for a sum separable potential is very similar but more complicated to analyse. 
The initial background trajectory must be similarly fine-tuned to be nearly parallel to one of the axes of the inflaton 



15 



0.15 0.5 0.85 1.0 




sinV 
0.15 0.5 0.85 1.0 



.sinV 



0.15 0.5 0.85 




0.85 




0.15 



-10 



-5 



FIG. 6: The function F as defined in Eqn. (JT9J) is plotted for three values of X: (a) X = -0.001 (left), (b) X = (middle) and 
(c) X = 0.001 (right). The central graph corresponds to f p since X — 0. 



fields, but the tuning at the end of inflation is not written so simply. We have also shown how the constraints on 
generating a large non-Gaussianity may be eased if we generalise the inflaton model to allow a non-canonical kinetic 
term. 

We have presented two explicit models where a large non-Gaussianity can be generated during slow roll inflation. 
One is a product potential driven by a field with a quadratic potential, which ends inflation when it approaches the 
minimum of the potential. If the second field with an exponential potential has a sufficiently small initial value then 
for certain values of the model parameter, A, this field grows by the right amount to generate a large non-Gaussianity, 
/nl ~ 100. We also consider the sum separable model of hybrid inflation. We provide a few specific choices of the 
model parameters, which are effectively r\ vlp and rj xx , for a suitably large ratio of the initial field values this model 
generates a large non-Gaussianity. We find this is possible either if the r/'s have the opposite sign (so that inflation 
takes place near a saddle point) or when both of the rfs have the same sign. However we have also found that if one of 
the field values is too small then quantum fluctuations may perturb the background trajectory to an extent that the 
classical trajectory in field space is no longer valid. There is also the possibility that the large scale loop correction 
which we have not generally considered in this paper may not be negligible. We intend to return to these issues and 
make a more thorough investigation of the hybrid model in a future publication. 
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